function [Tdt uxdt uydt P1dt P2dt] = metric(T, ux, uy, P1, P2, Tdx, Tdy, uxdx, uxdy, uydx, uydy, P1dx, P1dy, P2dx, P2dy)

%   choice of Pi: 
%   Pi_xx =  2 P1 (1+UX^2) +  2 UX UY P2
%   Pi_yy = -2 P1 (1+UY^2) +  2 UX UY P2
%   Pi_xy =  P2 ( 2 + UX^2 + UY^2 )
    global A;

    m11 = (-3).*T.^2.*(1+ux.^2+uy.^2).^(1/2);
    m12 = (1/2).*(1+ux.^2+uy.^2).^(-1/2).*((-3).*T.^3.*ux+(-4).*P1.*ux.*(1+uy.^2)+(-2).*P2.*uy.*(2+(-1).*ux.^2+uy.^2));
    m13 = (1/2).*(1+ux.^2+uy.^2).^(-1/2).*(((-3).*T.^3+4.*P1.*(1+ux.^2)).*uy+(-2).*P2.*ux.*(2+ux.^2+(-1).*uy.^2));

    m22 = (1/4).*(1+ux.^2+uy.^2).^(-3/2).*(2.*P2.*ux.*uy.*(2+3.*ux.^2+3.*ux.^4+uy.^2+(-3).*uy.^4)+3.*T.^3.*(ux.^4+3.*ux.^2.*(1+uy.^2)+2.*(1+uy.^2).^2)+4.*P1.*(ux.^4.*(1+(-3).*uy.^2)+2.*(1+uy.^2)+(-3).*ux.^2.*((-1)+uy.^4)));
    m23 = (1/4).*(1+ux.^2+uy.^2).^(-3/2).*(ux.*uy.*((-3).*T.^3.*(1+ux.^2+uy.^2)+4.*P1.*(1+ux.^2).*(1+3.*ux.^2+3.*uy.^2))+P2.*(8+(-6).*ux.^4+(-6).*ux.^6+12.*uy.^2+8.*uy.^4+ux.^2.*(8+6.*uy.^2+6.*uy.^4)));
    m24 = 2.*(ux+ux.^3).*(1+ux.^2+uy.^2).^(-1/2);
    m25 = uy.*(1+ux.^2+uy.^2).^(-1/2).*(2+3.*ux.^2+uy.^2);


    m32 = (1/4).*(1+ux.^2+uy.^2).^(-3/2).*(2.*P2.*(4+4.*uy.^2+(-3).*uy.^4+(-3).*uy.^6+3.*ux.^2.*(2+uy.^2)+ux.^4.*(4+3.*uy.^2))+(-1).*ux.*uy.*(3.*T.^3.*(1+ux.^2+uy.^2)+4.*P1.*(1+uy.^2).*(1+3.*ux.^2+3.*uy.^2)));
    m33 = (1/4).*(1+ux.^2+uy.^2).^(-3/2).*(2.*P2.*ux.*uy.*(2+ux.^2+(-3).*ux.^4+3.*uy.^2+3.*uy.^4)+3.*T.^3.*(2+2.*ux.^4+3.*uy.^2+uy.^4+ux.^2.*(4+3.*uy.^2))+4.*P1.*((-2)+(-3).*uy.^2+3.*ux.^4.*uy.^2+(-1).*uy.^4+ux.^2.*((-2)+3.*uy.^4)));
    m34 = (-2).*(1+ux.^2+uy.^2).^(-1/2).*(uy+uy.^3);
    m35 = ux.*(1+ux.^2+uy.^2).^(-1/2).*(2+ux.^2+3.*uy.^2);

    m42 = (1/12).*T.^(-4).*(1+ux.^2+uy.^2).^(-1/2).*(16.*(1+A).*P1.^2.*ux.*(1+uy.^2).*(2+ux.^2+uy.^2)+3.*T.^3.*(T.^3.*ux.*(2+ux.^2+3.*uy.^2)+2.*P2.*uy.*((2+A).*ux.^4+(-2).*ux.^2.*((-1)+A+2.*uy.^2+A.*uy.^2)+(2+uy.^2).*(2+(2+A).*uy.^2)))+2.*P1.*(4.*(1+A).*P2.*uy.*((-1).*ux.^4+(2+uy.^2).^2)+3.*T.^3.*ux.*(6+7.*uy.^2+4.*uy.^4+ux.^2.*(3+(-4).*uy.^2)+A.*(6+7.*uy.^2+2.*uy.^4+ux.^2.*(3+(-2).*uy.^2)))));
    m43 = (1/12).*T.^(-4).*(1+ux.^2+uy.^2).^(-1/2).*((-16).*(1+A).*P1.^2.*(1+ux.^2).*uy.*(2+ux.^2+uy.^2)+(-3).*T.^3.*(T.^3.*uy.*(2+3.*ux.^2+uy.^2)+2.*P2.*ux.*(4+(2+A).*ux.^4+(-2).*((-1)+A).*uy.^2+(2+A).*uy.^4+ux.^2.*(6+(-4).*uy.^2+(-2).*A.*((-1)+uy.^2))))+2.*P1.*(4.*(1+A).*P2.*ux.*(4+4.*ux.^2+ux.^4+(-1).*uy.^4)+3.*T.^3.*uy.*(4.*ux.^4+ux.^2.*(7+(-4).*uy.^2)+3.*(2+uy.^2)+A.*(2.*ux.^4+ux.^2.*(7+(-2).*uy.^2)+3.*(2+uy.^2)))));
    m44 = (1+A).*T.^(-1).*(1+ux.^2+uy.^2).^(1/2).*(2+ux.^2+uy.^2);


    m52 = (1/6).*T.^(-4).*(1+ux.^2+uy.^2).^(-1/2).*(3.*T.^6.*uy.*(1+uy.^2)+4.*(1+A).*P2.^2.*uy.*((-1).*ux.^4+(2+uy.^2).^2)+4.*P1.*(1+uy.^2).*((-3).*T.^3.*(1+(2+A).*ux.^2).*uy+2.*(1+A).*P2.*ux.*(2+ux.^2+uy.^2))+3.*P2.*T.^3.*ux.*(6+uy.^2+(-4).*uy.^4+ux.^2.*(5+4.*uy.^2)+A.*(6+uy.^2+(-2).*uy.^4+ux.^2.*(5+2.*uy.^2))));
    m53 = (1/6).*T.^(-4).*(1+ux.^2+uy.^2).^(-1/2).*(4.*(1+A).*P2.^2.*ux.*(4+4.*ux.^2+ux.^4+(-1).*uy.^4)+3.*T.^3.*ux.*(1+ux.^2).*(T.^3+4.*P1.*(1+(2+A).*uy.^2))+(-1).*P2.*uy.*(8.*(1+A).*P1.*(1+ux.^2).*(2+ux.^2+uy.^2)+3.*T.^3.*((-6)+4.*ux.^4+(-5).*uy.^2+(-1).*ux.^2.*(1+4.*uy.^2)+A.*((-6)+2.*ux.^4+(-5).*uy.^2+(-1).*ux.^2.*(1+2.*uy.^2)))));
    m55 = (1+A).*T.^(-1).*(1+ux.^2+uy.^2).^(1/2).*(2+ux.^2+uy.^2);

    
    R1 = (-3).*T.^2.*Tdx.*ux+(-3).*T.^2.*Tdy.*uy+(1/2).*uxdx.*(1+ux.^2+uy.^2).^(-1).*(2.*P2.*ux.*(ux+(-1).*uy).*uy.*(ux+uy)+(-4).*P1.*(1+ux.^2).*(1+uy.^2)+(-3).*T.^3.*(1+ux.^2+uy.^2))+(-1).*uxdy.*(1+ux.^2+uy.^2).^(-1).*(2.*P1.*ux.*uy.*(1+uy.^2)+P2.*(2+ux.^2+(-1).*((-3)+ux.^2).*uy.^2+uy.^4))+(1+ux.^2+uy.^2).^(-1).*(2.*P1.*ux.*(1+ux.^2).*uy+(-1).*P2.*(2+ux.^4+uy.^2+(-1).*ux.^2.*((-3)+uy.^2))).*uydx+(1/2).*(1+ux.^2+uy.^2).^(-1).*(4.*P1.*(1+ux.^2).*(1+uy.^2)+2.*P2.*ux.*uy.*((-1).*ux.^2+uy.^2)+(-3).*T.^3.*(1+ux.^2+uy.^2)).*uydy;
    R2 = (3/2).*T.^2.*Tdx+2.*P1dx.*(1+ux.^2)+2.*P2dx.*ux.*uy+P2dy.*(2+ux.^2+uy.^2)+(1/4).*uxdx.*(1+ux.^2+uy.^2).^(-1).*(3.*T.^3.*ux.*(1+ux.^2+uy.^2)+4.*P1.*ux.*(1+uy.^2+ux.^2.*(1+(-3).*uy.^2))+2.*P2.*uy.*(ux.^2.*(4+3.*ux.^2+(-3).*uy.^2)+4.*(1+uy.^2)))+(1/2).*uxdy.*(1+ux.^2+uy.^2).^(-1).*(P2.*ux.*((-2)+ux.^2+((-5)+3.*ux.^2).*uy.^2+(-3).*uy.^4)+3.*uy.*((-2).*P1.*ux.^2.*(1+uy.^2)+T.^3.*(1+ux.^2+uy.^2)))+(1/2).*ux.*(1+ux.^2+uy.^2).^(-1).*(6.*P1.*ux.*(1+ux.^2).*uy+P2.*((-2)+(-3).*ux.^4+uy.^2+ux.^2.*((-5)+3.*uy.^2))).*uydx+(1/4).*(1+ux.^2+uy.^2).^(-1).*(12.*P1.*(ux+ux.^3).*(1+uy.^2)+(-3).*T.^3.*ux.*(1+ux.^2+uy.^2)+2.*P2.*uy.*(4+4.*ux.^2+(-3).*ux.^4+(4+3.*ux.^2).*uy.^2)).*uydy;

    R3 = (3/2).*T.^2.*Tdy+2.*P2dy.*ux.*uy+(-2).*P1dy.*(1+uy.^2)+P2dx.*(2+ux.^2+uy.^2)+(1/2).*uxdy.*uy.*(1+ux.^2+uy.^2).^(-1).*((-6).*P1.*ux.*uy.*(1+uy.^2)+P2.*((-2)+ux.^2+((-5)+3.*ux.^2).*uy.^2+(-3).*uy.^4))+(1/4).*uxdx.*(1+ux.^2+uy.^2).^(-1).*((-3).*uy.*(4.*P1.*(1+ux.^2).*(1+uy.^2)+T.^3.*(1+ux.^2+uy.^2))+2.*P2.*ux.*(4+4.*uy.^2+(-3).*uy.^4+ux.^2.*(4+3.*uy.^2)))+(1/2).*(1+ux.^2+uy.^2).^(-1).*(3.*T.^3.*ux.*(1+ux.^2+uy.^2)+uy.*(6.*P1.*ux.*(1+ux.^2).*uy+P2.*((-2)+uy.^2+ux.^2.*((-5)+(-3).*ux.^2+3.*uy.^2)))).*uydx+(1/4).*(1+ux.^2+uy.^2).^(-1).*((-4).*P1.*(1+ux.^2).*uy+4.*P1.*((-1)+3.*ux.^2).*uy.^3+3.*T.^3.*uy.*(1+ux.^2+uy.^2)+2.*P2.*ux.*(4+4.*uy.^2+3.*uy.^4+ux.^2.*(4+(-3).*uy.^2))).*uydy;

    R4 = P1.*(2+ux.^2+uy.^2)+(1+A).*P1dx.*T.^(-1).*ux.*(2+ux.^2+uy.^2)+(1+A).*P1dy.*T.^(-1).*uy.*(2+ux.^2+uy.^2)+(1/6).*T.^(-4).*uxdy.*(1+ux.^2+uy.^2).^(-1).*(8.*(1+A).*P1.^2.*ux.*uy.*(1+uy.^2).*(2+ux.^2+uy.^2)+3.*T.^3.*((-1).*A.*P2.*(2+ux.^2).^2+T.^3.*ux.*(1+ux.^2).*uy+P2.*(4+(-4).*A+2.*ux.^2+(2+A).*ux.^4).*uy.^2+T.^3.*ux.*uy.^3+P2.*(6+A+(-2).*(2+A).*ux.^2).*uy.^4+(2+A).*P2.*uy.^6)+2.*P1.*(3.*T.^3.*uy.*(A.*ux.^3+(-1).*(2+A).*ux.*((-1)+ux.^2).*uy.^2+(2+A).*ux.*uy.^4)+2.*(1+A).*P2.*((2+ux.^2).^2+(8+2.*ux.^2+(-1).*ux.^4).*uy.^2+5.*uy.^4+uy.^6)))+(1/12).*T.^(-4).*uxdx.*(1+ux.^2+uy.^2).^(-1).*(16.*(1+A).*P1.^2.*(1+ux.^2).*(1+uy.^2).*(2+ux.^2+uy.^2)+2.*P1.*((-4).*(1+A).*P2.*ux.*(ux+(-1).*uy).*uy.*(ux+uy).*(2+ux.^2+uy.^2)+3.*T.^3.*(3.*(1+A).*(2+3.*ux.^2+ux.^4)+(9.*(1+A)+10.*(1+A).*ux.^2+(-2).*(2+A).*ux.^4).*uy.^2+(3.*(1+A)+2.*(2+A).*ux.^2).*uy.^4))+3.*T.^3.*(T.^3.*(1+ux.^2+uy.^2).*(2+ux.^2+uy.^2)+2.*P2.*ux.*uy.*((2+A).*ux.^4+(2+uy.^2).*(2+(2+A).*uy.^2)+(-2).*ux.^2.*((-1)+A+(2+A).*uy.^2))))+(1/6).*T.^(-4).*(1+ux.^2+uy.^2).^(-1).*((-8).*(1+A).*P1.^2.*ux.*(1+ux.^2).*uy.*(2+ux.^2+uy.^2)+2.*P1.*(2.*(1+A).*P2.*(2+ux.^2+uy.^2).*(2+ux.^4+uy.^2+(-1).*ux.^2.*((-3)+uy.^2))+3.*T.^3.*ux.*uy.*((2+A).*ux.^4+A.*uy.^2+(-1).*(2+A).*ux.^2.*((-1)+uy.^2)))+(-3).*T.^3.*(A.*P2.*((-2)+ux.^2+(-1).*uy.^2).*(2+ux.^4+uy.^2+(-1).*ux.^2.*((-3)+uy.^2))+ux.*(T.^3.*uy.*(1+ux.^2+uy.^2)+2.*P2.*ux.*(2+ux.^4+uy.^2+uy.^4+ux.^2.*(3+(-2).*uy.^2))))).*uydx+(1/12).*T.^(-4).*(1+ux.^2+uy.^2).^(-1).*((-16).*(1+A).*P1.^2.*(1+ux.^2).*(1+uy.^2).*(2+ux.^2+uy.^2)+(-3).*T.^3.*(T.^3.*(1+ux.^2+uy.^2).*(2+ux.^2+uy.^2)+2.*P2.*ux.*uy.*(4+2.*(3+A).*ux.^2+(2+A).*ux.^4+(-2).*((-1)+A+(2+A).*ux.^2).*uy.^2+(2+A).*uy.^4))+2.*P1.*(4.*(1+A).*P2.*ux.*(ux+(-1).*uy).*uy.*(ux+uy).*(2+ux.^2+uy.^2)+3.*T.^3.*(3.*(1+A).*(2+3.*ux.^2+ux.^4)+(9.*(1+A)+10.*(1+A).*ux.^2+2.*(2+A).*ux.^4).*uy.^2+(3.*(1+A)+(-2).*(2+A).*ux.^2).*uy.^4))).*uydy;
    R5 = P2.*(2+ux.^2+uy.^2)+(1+A).*P2dx.*T.^(-1).*ux.*(2+ux.^2+uy.^2)+(1+A).*P2dy.*T.^(-1).*uy.*(2+ux.^2+uy.^2)+(1/6).*T.^(-4).*uxdy.*(1+ux.^2+uy.^2).^(-1).*(3.*T.^3.*(1+uy.^2).*((4.*A.*P1+T.^3).*(1+ux.^2)+((-4).*P1+T.^3+(-4).*(2+A).*P1.*ux.^2).*uy.^2)+(-4).*(1+A).*P2.^2.*(2+ux.^2+uy.^2).*(ux.^2.*((-1)+uy.^2)+(-1).*(1+uy.^2).*(2+uy.^2))+2.*P2.*ux.*uy.*(4.*(1+A).*P1.*(1+uy.^2).*(2+ux.^2+uy.^2)+3.*T.^3.*(ux+(-1).*uy).*(ux+uy).*(1+(2+A).*uy.^2)))+(1/6).*T.^(-4).*uxdx.*(1+ux.^2+uy.^2).^(-1).*(4.*P1.*(1+uy.^2).*((-3).*T.^3.*ux.*(1+(2+A).*ux.^2).*uy+2.*(1+A).*P2.*(1+ux.^2).*(2+ux.^2+uy.^2))+P2.*((-4).*(1+A).*P2.*ux.*(ux+(-1).*uy).*uy.*(ux+uy).*(2+ux.^2+uy.^2)+3.*T.^3.*((1+A).*(6+9.*ux.^2+5.*ux.^4)+(9.*(1+A)+4.*(1+A).*ux.^2+2.*(2+A).*ux.^4).*uy.^2+(3.*(1+A)+(-2).*(2+A).*ux.^2).*uy.^4)))+(1/6).*T.^(-4).*(1+ux.^2+uy.^2).^(-1).*(4.*(1+A).*P2.^2.*(2+ux.^2+uy.^2).*(2+ux.^4+uy.^2+(-1).*ux.^2.*((-3)+uy.^2))+(-2).*P2.*ux.*uy.*(3.*T.^3.*(1+(2+A).*ux.^2).*(ux+(-1).*uy).*(ux+uy)+4.*(1+A).*P1.*(1+ux.^2).*(2+ux.^2+uy.^2))+3.*T.^3.*(1+ux.^2).*(T.^3.*(1+ux.^2+uy.^2)+4.*P1.*ux.^2.*(1+2.*uy.^2)+4.*A.*P1.*((-1)+((-1)+ux.^2).*uy.^2))).*uydx+(1/6).*T.^(-4).*(1+ux.^2+uy.^2).^(-1).*(4.*P1.*(1+ux.^2).*((-2).*(1+A).*P2.*(1+uy.^2).*(2+ux.^2+uy.^2)+3.*T.^3.*ux.*uy.*(1+(2+A).*uy.^2))+P2.*(4.*(1+A).*P2.*ux.*(ux+(-1).*uy).*uy.*(ux+uy).*(2+ux.^2+uy.^2)+3.*T.^3.*(3.*(1+A).*(2+3.*ux.^2+ux.^4)+(9.*(1+A)+4.*(1+A).*ux.^2+(-2).*(2+A).*ux.^4).*uy.^2+(5.*(1+A)+2.*(2+A).*ux.^2).*uy.^4))).*uydy;

    mdet = m11.*m25.*m34.*m43.*m52+(-1).*m11.*m24.*m35.*m43.*m52+(-1).*m11.*m25.*m33.*m44.*m52+m11.*m23.*m35.*m44.*m52+(-1).*m11.*m25.*m34.*m42.*m53+m11.*m24.*m35.*m42.*m53+m11.*m25.*m32.*m44.*m53+(-1).*m11.*m22.*m35.*m44.*m53+(-1).*m11.*m24.*m33.*m42.*m55+m11.*m23.*m34.*m42.*m55+m11.*m24.*m32.*m43.*m55+(-1).*m11.*m22.*m34.*m43.*m55+(-1).*m11.*m23.*m32.*m44.*m55+m11.*m22.*m33.*m44.*m55;

   
    i11 = m35.*((-1).*m24.*m43.*m52+m23.*m44.*m52+m24.*m42.*m53+(-1).*m22.*m44.*m53)+m25.*(m34.*m43.*m52+(-1).*m33.*m44.*m52+(-1).*m34.*m42.*m53+m32.*m44.*m53)+((-1).*m24.*m33.*m42+m23.*m34.*m42+m24.*m32.*m43+(-1).*m22.*m34.*m43+(-1).*m23.*m32.*m44+m22.*m33.*m44).*m55;
    i12 = (-1).*m13.*(m35.*m44.*m52+m34.*m42.*m55+(-1).*m32.*m44.*m55)+m12.*(m35.*m44.*m53+m34.*m43.*m55+(-1).*m33.*m44.*m55);
    i13 = m25.*m44.*(m13.*m52+(-1).*m12.*m53)+(m13.*m24.*m42+(-1).*m12.*m24.*m43+(-1).*m13.*m22.*m44+m12.*m23.*m44).*m55;
    i14 = (-1).*(m25.*m34+(-1).*m24.*m35).*(m13.*m52+(-1).*m12.*m53)+((-1).*m13.*m24.*m32+m12.*m24.*m33+m13.*m22.*m34+(-1).*m12.*m23.*m34).*m55;
    i15 = (m25.*m34+(-1).*m24.*m35).*(m13.*m42+(-1).*m12.*m43)+((-1).*m13.*m25.*m32+m12.*m25.*m33+m13.*m22.*m35+(-1).*m12.*m23.*m35).*m44;
    
    i22 = (-1).*m11.*(m35.*m44.*m53+m34.*m43.*m55+(-1).*m33.*m44.*m55);
    i23 = m11.*(m25.*m44.*m53+m24.*m43.*m55+(-1).*m23.*m44.*m55);
    i24 = m11.*((-1).*m25.*m34.*m53+m24.*m35.*m53+(-1).*m24.*m33.*m55+m23.*m34.*m55);
    i25 = m11.*(m25.*m34.*m43+(-1).*m24.*m35.*m43+(-1).*m25.*m33.*m44+m23.*m35.*m44);
    
    i32 = m11.*(m35.*m44.*m52+m34.*m42.*m55+(-1).*m32.*m44.*m55);
    i33 = (-1).*m11.*(m25.*m44.*m52+m24.*m42.*m55+(-1).*m22.*m44.*m55);
    i34 = m11.*(m25.*m34.*m52+(-1).*m24.*m35.*m52+m24.*m32.*m55+(-1).*m22.*m34.*m55);
    i35 = m11.*((-1).*m25.*m34.*m42+m24.*m35.*m42+m25.*m32.*m44+(-1).*m22.*m35.*m44);

    i42 = m11.*((-1).*m35.*m43.*m52+m35.*m42.*m53+(-1).*m33.*m42.*m55+m32.*m43.*m55);
    i43 = m11.*(m25.*m43.*m52+(-1).*m25.*m42.*m53+m23.*m42.*m55+(-1).*m22.*m43.*m55);
    i44 = m11.*((-1).*m25.*m33.*m52+m23.*m35.*m52+m25.*m32.*m53+(-1).*m22.*m35.*m53+(-1).*m23.*m32.*m55+m22.*m33.*m55);
    i45 = m11.*(m25.*m33.*m42+(-1).*m23.*m35.*m42+(-1).*m25.*m32.*m43+m22.*m35.*m43);

    i52 = m11.*(m34.*m43.*m52+(-1).*m33.*m44.*m52+(-1).*m34.*m42.*m53+m32.*m44.*m53);
    i53 = m11.*((-1).*m24.*m43.*m52+m23.*m44.*m52+m24.*m42.*m53+(-1).*m22.*m44.*m53);
    i54 = m11.*(m24.*m33.*m52+(-1).*m23.*m34.*m52+(-1).*m24.*m32.*m53+m22.*m34.*m53);
    i55 = m11.*((-1).*m24.*m33.*m42+m23.*m34.*m42+m24.*m32.*m43+(-1).*m22.*m34.*m43+(-1).*m23.*m32.*m44+m22.*m33.*m44);

    i11 = i11./mdet;
    i12 = i12./mdet;
    i13 = i13./mdet;
    i14 = i14./mdet;
    i15 = i15./mdet;

    i22 = i22./mdet;
    i23 = i23./mdet;
    i24 = i24./mdet;
    i25 = i25./mdet;

    i32 = i32./mdet;
    i33 = i33./mdet;
    i34 = i34./mdet;
    i35 = i35./mdet;

    i42 = i42./mdet;
    i43 = i43./mdet;
    i44 = i44./mdet;
    i45 = i45./mdet;

    i52 = i52./mdet;
    i53 = i53./mdet;
    i54 = i54./mdet;
    i55 = i55./mdet;

    
    Tdt  = -(i11.*R1+i12.*R2+i13.*R3+i14.*R4+i15.*R5);
    uxdt = -(i22.*R2+i23.*R3+i24.*R4+i25.*R5);
    uydt = -(i32.*R2+i33.*R3+i34.*R4+i35.*R5);
    P1dt = -(i42.*R2+i43.*R3+i44.*R4+i45.*R5);
    P2dt = -(i52.*R2+i53.*R3+i54.*R4+i55.*R5);




end
